Ground-state and dynamical properties of two-dimensional dipolar Fermi liquids 
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We study the ground-state properties of a two-dimensional spin-polarized fluid of dipolar fermions 
within the Eulcr-Lagrange Fermi-hypernetted-chain approximation. Our method is based on the 
solution of a scattering Schrodinger equation for the "pair amplitude" sjg(r), where g(r) is the pair 
distribution function. A key ingredient in our theory is the effective pair potential, which includes a 
bosonic term from Jastrow-Feenberg correlations and a fermionic contribution from kinetic energy 
and exchange, which is tailored to reproduce the Hartree-Fock limit at weak coupling. Very good 
agreement with recent results based on quantum Monte Carlo simulations is achieved over a wide 
range of coupling constants up to the liquid-to-crystal quantum phase transition (QPT). Using a 
certain approximate model for the dynamical density-density response function, we furthermore 
demonstrate that: i) the liquid phase is stable towards the formation of density waves up to the 
liquid-to-crystal QPT and ii) an undamped zero-sound mode exists for any value of the interaction 
strength, down to infinitesimally weak couplings. 

PACS numbers: 03.75.Ss, 67.85.-d, 67.85.Lm 



I. INTRODUCTION 

Recent experimental breakthroughs in trapping and 
cooling polar molecules and atoms with large permanent 
magnetic moments has triggered an immense theoretical 
interest in quantum dipolar fluids [IHS] . Unlike the usual 
van der Waals interaction between atoms, which can be 
replaced by a contact Fermi pseudo-potential at ultra- 
low temperatures [5] , the dipole-dipole interaction is long 
ranged and anisotropic. It is therefore natural to expect 
more exotic phases in ultra-cold dipolar gases. While 
one of the greatest advantages of short-range interactions 
is their tunability through Feshbach resonances [BJ [TJ, 
techniques have been proposed [Sj for controlling both 
strength and sign of dipolar interactions as well. 

As already mentioned, the inter-particle interaction be- 
tween polarized (i.e. dipoles aligned in the same direc- 
tion) dipoles has two important features: i) it is long- 
ranged, i.e. it decays like 1/r 3 at large distances, and ii) 
it is anisotropic. In particular, it is repulsive for dipoles 
aligned side-by-side and is attractive for dipoles aligned 
head-to-toe. 

It is worth mentioning that the attractive part of 
dipole-dipole interactions can drive dipolar fluids towards 
instabilities. In alkali- metal diatomic molecules such as 
K-Rb, Li-Na, etc., some chemically reactive channels, 
which are energetically favorable, exist and lead to par- 
ticle recombination and two-body losses in the gas [5] . 

A very simple method for stabilizing dipolar gases is 
to confine them into low-dimensional geometries. For 
example, a trap with pancake geometry together with 
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a polarizing field, which aligns the dipoles along the 
direction of transverse confinement, simulates a stable 
two-dimensional (2D) system with purely repulsive and 
isotropic dipolar interactions of the form 



Here Cdd is the dipole-dipole coupling constant, which 
depends on the microscopic origin of the interaction: e.g., 
it is d 2 /eo for particles with permanent electric dipole d 
and [IqM 2 for particles with permanent magnetic dipole 
M (here eo and /io are the permittivity and permeability 
of vacuum, respectively). 

Ground-state properties and collective modes of 2D 
dipolar fermions have been addressed in a number of 
studies (TUH22] • For their particular relevance to this Ar- 
ticle we highlight the following two recent studies [3TJ [55] 
of a 2D dipolar Fermi gas (DFG) with isotropic interac- 
tions as in Eq. Q. Lu and Shlyapnikov [21] have cal- 
culated a number of Fermi-liquid properties of a weakly 
interacting 2D DFG. In particular, these authors have 
presented several exact results up to second order in a 
natural dimensionless coupling constant, which we have 
introduced below in Eq. Q. More recently, Matveeva 
and Giorgini [55] have carried out quantum Monte Carlo 
(QMC) simulations of a 2D DFG, presenting in partic- 
ular results for the phase diagram of this system over 
a wide range of coupling constants. These studies pose 
severe bounds on any microscopic theory of 2D DFGs. 

In this Article we present a theoretical study of ground- 
state and dynamical properties of a 2D DFG with av- 
erage density n. Our main focus is on the pair dis- 
tribution function (PDF) g(r), which is often referred 
to as "Pauli-Coulomb hole". This is defined [23] [24] 
so that the quantity 2irng(r)dr gives the average num- 
ber of dipoles lying within a circular shell of radius 
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dr centered on a "reference" dipole sitting at the ori- 
gin. We present a self-consistent semi-analytic theory of 
the PDF, which incorporates many-body exchange and 
correlation effects, thereby allowing us to explore the 
physics of the system at strong coupling. Our approach, 
which is based on the so-called Euler-Lagrange Fermi- 
hypernetted-chain (FHNC) approximation at zero tem- 
perature |25H28j . involves the solution of a zero-energy 
scattering Schrodinger equation with a suitable effective 
potential [29Tl34j . This contains a "bosonic term" from 
Jastrow-Feenberg correlations and a "fermionic term" 
from kinetic energy and exchange, which is tailored to 
reproduce the Hartree-Fock (HF) limit at weak coupling 
and guarantees the antisymmetry of the fermionic wave 
function. Furthermore, we use the fluctuation-dissipation 
theorem [21 |2j|] and the PDF obtained from the FHNC 
approximation to calculate the dynamical density-density 
linear-response function. With this quantity at our dis- 
posal, we investigate the possibility of instabilities to- 
wards inhomogenous ground states (i.e. density waves) 
at strong coupling and the existence of a "zero sound" 
mode [23] in a 2D DFG. Our results are severely bench- 
marked against the findings of Refs. [2T1 122) . 

This Article is organized as follows. In Sect. [IT] we 
present our model and the self-consistent method we use 
to calculate in an accurate manner the PDF of a 2D 
DFG. In Sect. [TlT] we discuss a number of approximations 
we make to derive the dynamical density-density linear- 
response function of a 2D DFG and explain how this can 
be used to examine the tendency towards a density-wave 
instability and the emergence of a collective zero-sound 
mode due to many-body effects. Sect. IV collects our 
main numerical results, while Sect. [V] contains a brief 
summary of our main findings. 



II. SCATTERING THEORY FOR THE 
PAULI-COULOMB HOLE 

We consider a spin-polarized 2D DFG described by the 
following first-quantized Hamiltonian [35) : 



u 



2m 



(2) 



i<3 



where m is the mass of a dipole and the bare dipolc- 
dipole interaction has been introduced above in Eq. ([lj. 
The ground-state properties of the Hamiltonian ([2| are 
governed by a single dimensionless parameter: 



A = fcpr 



(3) 



where ro = mCdd/(47r/i 2 ) is a characteristic length scale 
and kp — v Ann is the Fermi wave number, n being the 
2D average density. 

In order to calculate the ground-state properties of the 
Hamiltonian |2|, we use the FHNC [25H27] approxima- 
tion at zero temperature. In what follows we first present 
our theory at the simplest level (which works well in the 



perturbative regime A <C 1) and then transcend it to ob- 
tain accurate results at strong coupling (A ^> 1). 

With the zero of energy taken at the chemical poten- 
tial, one can write a formally exact differential equation 
for the PDF [31 [31]: 



-V 2 + V eS (r) 



(4) 



We write the effective scattering potential V e g(r) as the 
sum of three contributions: 



V eS (r) = v dd (r) + W B (r) + W F (r) 



(5) 



Here v dd (r) is the bare repulsive dipole-dipole interaction 
in Eq. while the bosonic contribution to the scatter- 
ing potential, Wsir), is defined, at the level of the so 
called "FHNC/0" approximation, by the following equa- 
tion 1551: 



Ws(fc) I FHNC/0 



e{k) 
2n 



[2S(k) + 1] 



S(k) - 1 

m 



(6) 



In writing Eq. (pj we have introduced the Fourier trans- 
form (FT) WB(fc) of Wb(?*) according to 



FT[F(r)] = J d 2 r F(r) exp (ik ■ r) 



(7) 



Furthermore, e(fc) = h 2 k 2 /(2m) is the single-particle en- 
ergy and S(k) is the instantaneous or "static" structure 
factor [23], S(k) = l + n FT[g(r) - 1]. 

When A ~ 1 the simplest approximation for Wb(t) in 
Eq. ([6| is inadequate. Improvements on Eq. ([6| can be 
sought in two directions [37] , The FHNC/0 may be tran- 
scended by the inclusion of (i) low-order "elementary" 
diagrams and (ii) three-body Jastrow-Feenberg correla- 
tions. 

The contribution from three-body correlations to the 
bosonic potential is given by [37) : 



Wi 3 \k) 



1 



d 2 q 



An 2 J (2tt) 2 

[E{p)+E(q)}u 3 (q,P,k)} 



S(p)S(q)u 3 (q,p,k){v 3 (q,p,k) 

(8) 



In the previous equation, p = —(q+k), E(k) = e(k)/S(k) 
is the Bijl-Feynman excitation spectrum |24j . 

v 3 (q,p,k) = (h 2 /m)[k-px(p) + k-q X (q)+p-qx(q)} , (9) 
and 

(h 2 /2m) 



u 3 (q,p,k) 



E(k) + E(p) + E{q) 
x [k ■ px(k)x(p) + P ■ Qx(p)x(l) 
+ k ■ qx(k) X (q)} ■ (10) 



In Eqs. (|9|-([T0]) xO) = 1 - l/S(k). We have taken 
into account higher-order terms that are missed by the 
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FHNC/O approximation by assuming that they lead to 
corrections to the scattering potential V e g(r). Using the 
theory developed by Apaja et al. [37], we have supple- 
mented We(fc) | FHNC/o m (|6j) by the inclusion of the 

/o\ 

three-body potential Hg (k): 



WB(k)\ FWC/a3 = W B (k)\ FliNC/0 + a(\)wi 3) (k) . ( 
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If a(A) is set to unity, the r.h.s. of Eq. (Ill de- 



fines the so-called "FHNC/3" approximation. It has 
been shown [32r 34j that higher-order corrections beyond 
FHNC/3 can be effectively taken into account by intro- 
ducing a weighting function a(X) > 1. This approxima- 
tion has been termed [31] "FHNC/a3". A convenient 
analytical parametrization of the function a(A) for 2D 
dipolar fluids can be found in Ref. [33]. Using the nota- 
tion of this Article, it reads as follows: 



where £o = £f/2 = ft 2 fcp/(4m) = fi 2 A 2 /(4mr 2 ) is the 
ground-state energy of the non-interacting system, £p be- 
ing the Fermi energy, and <7 7 (r) is the PDF of an auxil- 
iary system with scaled dipole-dipole interactions of the 

form u^d ( r ) — jVdd(f) = 7Cdd/(47rr 3 ). In practice, the 
integration over 7 is carried out by integrating over the 
coupling constant A. 

In Sect. |IV| we present numerical results obtained 
only within our most elaborate approximation, i.e. the 
FHNC/a3 approximation. Nevertheless, for the sake of 
simplicity, all our numerical results for g(r), S(k), and 
£gs w iU be labeled by the acronym "FHNC" (rather than 
"FHNC/a3"). 



III. LINEAR-RESPONSE THEORY, 
DENSITY- WAVE INSTABILITIES, AND 
COLLECTIVE MODES 



a(X) = 1.88 + 3.26exp (-0.046 A i ib ) 



(12) 



The previous equation is valid all the way up to the crit- 
ical coupling (A ~ 25) for the liquid-to-crystal quantum 
phase transition [2"2] . 

We finally turn to describe the last term in Eq. ([5]), 
which is supposed to take care of the fermionic statistics 
of the problem. According to the original version of the 
FHNC theory [25H27J, the "Fermi potential" W F (r) has 
a very complicated form. Here we have decided to use 
a simple but effective recipe, which was first proposed 
by Kallio and Piilo [29] for the 3D electron liquid. In 
this approximate scheme W F (r) is given by the following 
expression: 



W F (r) 



ft 2 VrvWfr) 



lim Wb(t) , 

A^O 



(13) 



where gnF(r) is the well-known [231 [23] 2D HF PDF and 
Wb{t) is the bosonic potential defined above in Eq. Q 
(at the FHNC/0 level) or in Eq.JTl) (at the FHNC/a3 
level). The simple choice in Eq. (|13[) guarantees that the 
HF limit is recovered exactly in weak coupling A — s- 
limit. The Fermi potential (13) has been extensively in- 
vestigated for 3D [30] and 2D|32j electron liquids yield- 
ing results in excellent agreement with QMC simulation 
data. 

Equations and (13) form a closed set of equa- 

tions, which can be solved numerically in a self-consistent 
manner to the desired degree of accuracy. Practical 
recipes on how to solve this system of equations are dis- 
cussed in detail in Ref. [30] . 

Once the PDF has been calculated, the ground-state 
energy per particle of the system, eqs, can be eas- 
ily extracted by using the integration-over-the-coupling- 
constant algorithm [21] : 



£GS = £0 



dj / d r v dd (r)g 7 (r) , 



(14) 



The density-density linear-response function of a 
many-particle system can be generically written as fol- 
lows [23]: 

Xnn(k,Uj) = - — — ——=———, (15) 

l — ip{k,u})xo{k,u3) e(k,u) 

where ijj(k,uj) is a suitable dynamical effective 
potential — not to be confused with the FT of the 
effective potential V c s(r) which enters the zero-energy 
scattering Schrddinger equation Q — and xo{k,oj) is the 
well-known [231 [35] density-density response function of 
an ideal (i.e. non-interacting ) 2D Fermi gas. 

In the celebrated Random Phase Approximation 
(RPA) [231 [24], the effective potential ip(k,u) is bru- 
tally approximated with the FT of the bare inter-particle 
potential, i.e. Udd(V) in our case. It is very well 
known [23J [24] that the RPA neglects short-range ex- 
change and correlation effects and that it is intrinsically 
a weak-coupling theory. It is thus not expected to work 
well (in reduced spatial dimensions and) for values of the 
dimensionless coupling constant A > 1. One of the main 
drawbacks of the RPA is that it grossly overestimates the 
strength of the Pauli-Coulomb hole by predicting large 
and negative values for g(r) at short distances, thereby 
violating the fundamental request g(r) > 0. Moreover, 
in the context of dipolar Fermi gases, the RPA predicts 
that the long-wavelength collective excitation spectrum 
(zero-sound mode) is empathic to the short-range details 
of the bare interaction potential [14"] . i.e. the ultraviolet 
cut-off which is needed to regularize the FT of the bare 
dipole-dipole potential VddM- 

In the past sixty years or so, a wide body of litera- 
ture has been devoted to transcend the RPA, especially 
in the context of 2D electron liquids in semiconductor 
heterojunctions [23]- Following the seminal works by 
Hubbard [39] and Singwi, Tosi, Land, and Sjdlander [40] 
(STLS), one successful route has been based on the use of 
"local field factors" [231 131] (LFFs) . Here we will not use 
Hubbard or STLS LFFs. (For a successful employment of 
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the STLS approximation in the context of 2D DFGs see, 
for example, Ref. pj5].) In this Article we would like to 
construct a reliable approximation for the density-density 
response function Xnn(k, lu), which is based on the FHNC 
theory of the PDF outlined in Sect. [IT] 

We thus start from the well-known fluctuation- 
dissipation theorem (FDT) [53] , which relates the imagi- 
nary part of density-density response function Xnn {k, ui) 
to the instantaneous structure factor S(k). At zero tem- 
perature the FDT reads [53] 



S{k) 



d,UJ 3m [Xnn{k,U))\ 



(16) 



To make some progress, we neglect the frequency de- 
pendence of the effective potential ip(k,Lu) in Eq. (15) 



we replace the complex function %p(k,tu) by a real quan- 
tity, which we denote by the symbol ip(k). This approx- 
imation is often made in treating correlation effects in 
the electron liquid [53] HI] an d is certainly shared by the 
most elementary theories based on LFFs (Hubbard and 
STLS). In this case, one can view Eq. (16) as an integral 
equation for the unknown quantity 4>(k), assuming that 
the l.h.s. of Eq. (16), i.e. the static structure factor, is 
accurately known, e.g. from QMC simulations or micro- 
scopic theories such as the one outlined in Sect. [Tl] This 
fully numerical approach has been successfully employed 
in different contexts [35J |33] ■ The physical interpretation 
of f/>(/c) is clear: it represents the "best" average effec- 
tive potential [averaged over frequency, as from Eq. ( 16 1] 
which, by virtue of the FDT, makes the response of the 
system consistent with the local structure of the fluid 
around a reference dipole (the Pauli-Coulomb hole). 

In the spirit of making the problem at hand more 
amenable to a semi-analytical treatment, we also use 
the so-called "mean-spherical approximation" (MSA) for 



(MSA),, s 



2ne{k) 



(hio + tO+y-[e(k)/S HF (k)Y 



(17) 



where Shf(&) is the well-known 2D HF static structure 
factor [24j . This approximation allows us to perform the 
integration over lu in Eq. ( 16 ) analytically, yielding 



4>(k) 



MSA 



£(*) 
2n 



S 2 (k) S^ F (k) 



(18) 



For the static structure factor S(k) in the r.h.s. of 
Eq. (fTsj) we use the FHNC theory described above in 
Sect.W 

We can now use Eqs. ( 15 ) and ( 18 ) to address two 
important issues. 

First, we can carry out a linear-stability analysis of 
the liquid phase against density modulations. In this re- 
spect, a pole in the static density-density response func- 
tion Xnn{k c ,to = 0) at a finite wave vector k c signals an 
instability of the liquid state towards a density wave with 



period oc k c 1 . In practice, we need to find whether the 
following equation, 



e(k, w = 0) = 1 - tp(k)xo(k, oj = 0) = 



(19) 



admits a solution at a finite wave vector k c . We remind 
the reader that, in the static limit, Xo(k, lu) is purely real. 

Second, we can study the existence of a collective 
mode [53] in the density channel (zero sound 23J). This 
is the solution of the complex equation e(k,Lu) = or, 
equivalently, of the following two real equations: 



l-$(k)3te [xo(k ) co)}=0 
3to [xo{k,co)] = 



(20) 



The solution luzs — ^zs(k) of Eq. (20) corresponds to 



a self-sustained oscillation with a non-trivial dispersion 
relation and a finite velocity v%s = lim^o w zs (k) /k in 
the long- wavelength limit. The second equation means 
that the collective mode is undamped when it falls in the 
region of (fc, lu) space where particle-hole pairs are absent. 
This occurs when vzs > vp, vp = hkp/m being the Fermi 
velocity. When the collective mode enters the particle- 
hole continuum, Landau damping starts: the mode has 
sufficient energy to decay by emitting a particle-hole pair 
while, at the same time, conserving momentum. 

Before concluding this Section, we derive a formal ex- 
pression for the ZS velocity vzs, in terms of V>(0). (As 
we will see below in Sect. IV ^p(k) is regular and positive 
at k = 0.) In order to find vzs we use the following long- 
wavelength limit of the ideal response function [53] [33] : 



lim 3?e [xo(k, vpku)] = -N(0) 1 - 



1 



(21) 



where N(0) = mj (2irh 2 ) is the 2D density-of-states at the 
Fermi energy and v = Lu/(vpk) = constant. Note that, 
in Eq. (21), the ratio between lu = vpkv and k remains 



constant in the limit k — > 0, precisely as in the ZS mode 
[lim/ £ _j.o LUzs(k)/k = constant]. It is very important to 
observe that the asymptotic behavior (21 ) needed for the 



calculation of the ZS velocity is very different from the 
usual high-frequency limit imposed by the f-sum rule [53] : 



5fe [xo{k,Lu > Vpk)] = 



nk 2 



(22) 



Now, replacing Eq. (21) [and not Eq. (22 



we find the following formal expression for t 
in units of the Fermi velocity: 

vzs = _1 + A(0)^(0) 
VF • / l + 2iV(0)^(0) 



in Eq. (20), 



le ZS velocity 



(23) 



which is well defined if tp(0) > -piVfO)]- 1 . Note that 
the quantity on the r.h.s. of Eq. ( [23] ) is always larger 
than one. We therefore conclude that, within the ap- 
proximations we made to derive Eq. ( 18 1, a 2D DFG dis- 
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FIG. 1: (Color online) The pair distribution function g(r) of 
a 2D fluid of dipolar fermions is plotted as a function of the 
scaled distance k?r and for various values of the dimensionless 
coupling constant A. In this plot, lines label the results of the 
FHNC approximation (this work) while symbols label QMC 
results [32]. 



A) an undamped ZS mode, in agreement with Ref. [2"T] . 
As discussed at length in Ref. [21], this mode stems en- 
tirely from correlation effects and it is thus not describ- 
able within the HF approximation. However, the RPA, 
which is the minimal theory including correlations, is not 
enough in this respect since it yields a ZS mode with a 
velocity that depends on the short-range cut off of the 
bare dipole-dipole interaction P3]. A serious theory of 
the ZS mode thus requires the inclusion of correlation 
effects beyond RPA. The FHNC theory discussed in this 
Article is an example. 



IV. NUMERICAL RESULTS AND DISCUSSION 

In this Section we present our main numerical results. 

We begin by showing our results for the PDF g{r) and 
static structure factor S(k). In Figs. [I] and [2] we com- 
pare our results (lines) with the corresponding QMC data 
(symbols) [55]. The agreement between theory and nu- 
merical simulations is clearly excellent up to very large 
values of the dimensionless coupling constant A (A = 20). 
At these values of A, conventional theories such as RPA 
and STLS fail even qualitatively. Note that, according to 
the QMC study by Matveeva and Giorgini [22] , a liquid- 
to-crystal quantum phase transition is expected to occur 
at A ~ 25. This is clearly signaled by the amplitude 
of the first-neighbor peak in the static structure factor 
(see Fig. [2]), which increases with increasing A indicat- 
ing the build up of correlations in the liquid phase upon 
approaching crystalline order. 

The PDF shown in Fig. [l] can be used to calculate the 
ground-state energy by employing Eq. ( 14 ) . In Fig. [3] 



FIG. 2: (Color online) The instantaneous structure factor 
S(k) of a 2D fluid of dipolar fermions is plotted as a function 
of fc/fcp and for various values of A. In this plot, lines label the 
results of the FHNC approximation (this work) while symbols 
label QMC results [22] . 
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FIG. 3: (Color online) The ground-state energy of a 2D 
fluid of dipolar fermions [in units of the Hartree-Fock en- 
ergy £hf/eo = 1 + 128A/(457r)], is plotted as a function of 
A. Blue circles label the results of the FHNC approximation 
(this work) while red squares label the QMC results |22j . The 
solid line represents the parametrization formula in Eq. | |24[ ) 
with a = 1.5006, b = 1.0107, and c = 0. The dashed line 
represents the parametrization in Eq. ( 24 1 with a — ( see 
main text), b = 1.1017, and c = -0.0100. 



the same plot we compare our findings with the corre- 
sponding QMC results [22]. In passing, we note that our 
FHNC results for the ground-state energy (per particle) 
can be accurately parametrized by the following expres- 
sion: 



£gs(A) 



£0 



128 

45"7t' 



we report our results for the ground-state energy as ob- 
tained from the PDF calculated at the FHNC level. In 



In 1 



1 



'X + b A + c A 3 / 2 , 



,(24) 
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where a, b, and c are numerical constants. The sum of 
the first two terms in square brackets on the r.h.s. of 
Eq. (24) yields the HF approximation for the ground- 
state energy [2T]: £hf = £o[l + 128A/(45tt)]. The best 
fit of our FHNC data for the energy of the liquid phase 
up to A = 40 is obtained by using a and b as free fitting 
parameters and setting c = 0: we find a = 1.5006 and 
b = 1.0107. The result of this two-parameter fit is shown 
in Fig. [3] (solid line). 

Alternatively, the simple formula in Eq. (24) can be 



used to parametrize also the QMC data by Matveeva 
and Giorgini [35]. Since these data are believed to be 
essentially exact, we can fix the value of a by imposing 
that Eq. (24) reproduces exactly the results of second- 
order perturbation theory |21j . Straightforward algebraic 



manipulations on Eq. ( 24 ) yield the following expansion 
in powers of A for A 



£gs(A) = £ 



0: 



128 X A W 2 XX 

A H ln(a 2 A) 



(25) 



where ". . . " denotes higher-order terms. To the same 
order of perturbation theory, Lu and Shlyapnikov [2Tj 
find [Eq. (91) in their work]: 




FIG. 4: (Color online) Three important thermodynamic 
quantities: the chemical potential fi, pressure P, and inverse 
compressibility « _1 of a 2D fluid of dipolar fermions (in units 
of their non-interacting values) are plotted as functions of A. 
Lines label analytic results obtained from the parametrization 
formula ( |24[ | while symbols label numerical results obtained 
directly from the FHNC ground-state energy. 



£GS(A) = £ 



^A+^ln(C 3 A) 

457T 4 



(26) 



in the 



where £3 = 1-43 (we have taken the limit A 
expression for (3 given in Ref. [21]). Comparing Eq. (25) 
with Eq. (26) we conclude that a 



1.2. The 



parameters b and c can then be used to yield the best fit 
to the QMC data for the energy of the liquid phase up 
to A = 72 [22]: we find b = 1.1017 and c = -0.0100. The 
result of this two-parameter fit is also shown in Fig. [3] 
(dashed line). 

The difference between the total ground-state energy 
and the non-interacting contribution e defines the inter- 
action energy: £i n t(A) = £qs(A) — £o(A). Note that unlike 
the gellium model for electron gases [23], the Hartree 
contribution to the interaction energy does not vanish 
in our system of polarized DFGs [21]. Eq. (24 1 thus 



provides an extremely useful input for calculations of 
ground-state properties of inhomogenous 2D DFGs based 
on density functional theory (DFT) [24]. In DFT, in- 
deed, one needs to approximate the unknown interac- 
tion energy E- mt [n(r)], viewed as a functional of the local 
ground-state density n(r). In the local density approxi- 
mation (LDA) one can write 24 



E int [n(r)] L ^ A J 



d 2 r n(r)s int (X(r)) , (27) 



where A(r) is defined as in Eq. ^ with n replaced by 
the local density n(r). An example where the DFT- 
LDA approach could be very useful is a 2D DFG in the 
presence of an in-plane harmonic confinement potential 



From the knowledge of the ground-state energy (per 
particle) £qs we can also construct a number of ther- 
modynamic quantities at zero temperature. Most no- 
tably, the chemical potential /1 = d(nscs)/dn, the pres- 
sure P = n 2 deGs/dn, and the inverse compressibility 
kT 1 = ndP/dn are readily obtained from the interpola- 



tion formula given in Eq. ( 24 ) . We display these quanti- 
ties as functions of the interaction strength A in Fig. [4] 
Note that all these quantities, which still remain to be 
experimentally measured, are strongly enhanced by in- 
teractions. 

Fig. [5] illustrates the effective potential i/>(fc) as ob- 
tained from Eq. (18). We clearly see from this plot that 
ip(k) is regular and positive at k = 0. 

In Fig. [6] we plot e{k, 0) = 1 — ^(k)xo(k,uj = 0) as a 
function of wave vector k and for different values of A. 
Increasing the interaction strength, a minimum occurs in 
e(fc,0) (yielding a peak in the density-density response 
function) at a wave vector close to 2fcp. This minimum 
remains finite, though, up to the largest value of A we 
have investigated (A = 40). In other words, our the- 
ory does not predict any density- wave instability in a 
2D DFG. This is in agreement with the QMC results by 
Matveeva and Giorgini [22] , who have shown that a stripe 
phase has higher energy than that of liquid and crystal 
phases at any A. 

Finally, in Fig. [7] we illustrate our predictions for the 
dispersion of the ZS mode. As already discussed at the 
end of Sect. |HI| our theory predicts an undamped ZS 
mode at long wavelengths for every value of A. The ZS 
velocity as well as the critical wave vector at which Lan- 
dau damping starts increase with increasing A. 



0.5 1 1.5 2 2.5 3 3.5 4 

k/k-p 



0.2 0.4 0.6 0.8 1 1.2 

fc/fcp 



FIG. 5: (Color online) The effective interaction tp(k) [in units 
of 2-kH 2 /m] in a 2D fluid of dipolar fermions as obtained from 
Eq. ( 18 1 is plotted as a function of k/k-p for various values of 
A. Note that i/)(fe — > 0) tends to a positive value. 
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FIG. 6: (Color online) The static dielectric function e(k, oj = 
0) of a 2D fluid of dipolar fermions as a function of k/kp and 
for various values of A. Note that e(k,oj = 0) never vanishes, 
even for very large values of A. 



V. SUMMARY 

In summary, we have presented an extensive study of 
ground-state and dynamical properties of a strongly cor- 
related two-dimensional spin-polarized fluid of dipolar 
fermions. 

The main focus of our work has been on the pair dis- 
tribution function g(r), a key ground-state property of 
any quantum fluid. To calculate the pair distribution 
function we have employed the Fermi-hypernetted-chain 
approximation combined with a zero-energy scattering 
Schrodinger equation for the "pair amplitude" yj g(f). 
The effective potential that enters this equation includes 
a bosonic term from Jastrow-Feenberg correlations and a 
fermionic contribution from kinetic energy and exchange, 



FIG. 7: (Color online) The frequency u>zs(k) of the zero-sound 
mode in a 2D fluid of dipolar fermions is plotted as a function 
of k/kp for various values of A. The thin (black) solid line rep- 
resents the upper bound of the particle-hole (p-h) continuum, 
i.e. cj+(k) = e(k)/H + vpk. Note that the zero-sound mode 
lies above the p-h continuum for every A, up to a A-dependent 
critical wave vector at which Landau damping starts. 



which is tailored to reproduce the Hartree-Fock limit 
at weak coupling. Our results for the pair distribution 
function and the static structure factor S(k) have been 
severely benchmarked against state-of-the-art quantum 
Monte Carlo results by Matveeva and Giorgini [22] . Very 
good agreement with these results has been achieved over 
a wide range of coupling constants up to the liquid-to- 
crystal quantum phase transition. 

By combining our knowledge on the pair distribu- 
tion function with the fluctuation-dissipation theorem, 
we have been able to calculate in an approximate fash- 
ion also the dynamical density-density response function. 
This ingredient has been used to demonstrate that, in a 
two-dimensional spin-polarized fluid of dipolar fermions, 
i) the liquid phase is stable towards the formation of 
density waves up to the liquid-to-crystal quantum phase 
transition (in agreement with Ref. [22]) and ii) an un- 
damped zero-sound mode occurs for any value of the 
interaction strength, down to infinitesimally weak cou- 
plings (in agreement with Ref. [21]). 

Last but not least, we have presented a useful 
parametrization formula, Eq. (24), for the ground-state 



energy of a two-dimensional spin-polarized fluid of dipo- 
lar fermions, which fits well both our Fermi-hypernetted- 
chain results and the quantum Monte Carlo data by 
Matveeva and Giorgini [22J. This can be very effectively 
employed in density-functional calculations of 2D inho- 
mogenous dipolar fermions. 
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